Symmetry Breaking in Symmetric and Asymmetric Double- Well Potentials 
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Motivated by recent experimental studies of matter-waves and optical beams in double well po- 
tentials, we study the solutions of the nonlinear Schrodinger equation in such a context. Using a 
Galerkin-type approach, we obtain a detailed handle on the nonlinear solution branches of the prob- 
lem, starting from the corresponding linear ones and predict the relevant bifurcations of solutions 
for both attractive and repulsive nonlinearities. The results illustrate the nontrivial differences that 
arise between the steady states/bifurcations emerging in symmetric and asymmetric double wells. 



Introduction. It is well known that the nonlinear 
Schrodinger (NLS) equation is a fundamental model de- 
scribing the evolution of a complex field envelope in non- 
linear dispersive media jlj. As such, it plays a key role 
in many different contexts, ranging from nonlinear and 
atom optics to plasma physics, fluid dynamics, and even 
biophysical models ^. The interest in the NLS equa- 
tion has dramatically increased during the last few years, 
as it also describes the mean-field dynamics of Bose- 
Einstein condensates (BECs) |3]- In this context, the 
NLS is also known as the Gross-Pitaevskii (GP) equa- 
tion, and typically incorporates external potentials that 
are used for the BEG confinement. Such potentials may 
be, e.g., harmonic (usually implemented by external mag- 
netic fields) or periodic (implemented by the interference 
of laser beams), so-called optical lattices |j]- Importantly, 
NLS models with similar external potentials appear also 
in the context of optics, where they respectively describe 
the evolution of an optical beam in a graded-index waveg- 
uide or in periodic waveguide arrays [^ @ . 

Another type of external potential, which has mainly 
been studied theoretically in the BEG context |7l.lMl3.[lOl 
[ll|,[l2| is the double well potential. Moreover, it has been 
demonstrated experimentally that a BEG either tunnels 
and performs Josephson oscillations between the wells, 
or is subject to macroscopic quantum self-trapping [l3l| . 
On the other hand, in the context of optics, a double 
well potential can be created by a two-hump self-guided 
laser beam in Kerr media [lj| . A different alternative was 
offered in jl^ , wherein the first stages of the evolution of 
an optical beam, initially focused between the wells of a 
photorefractive crystal, were monitored. 

One of the particularly interesting features of either 
matter-waves or optical beams in double well potentials 
is the spontaneous symmetry breaking, i.e., the localiza- 
tion of the respective order parameter in one of the wells 
of the potential. Symmetry breaking solutions of the NLS 
model have first been predicted in the context of molec- 
ular states |la | and, apart from the physical contexts 
of BEGs 0,11 H Ea O 113 (see also ll\) and optics 
[Tj, [i3 mentioned above, they have also been studied 
from a mathematical point of view in Refs. [l§.ll9|. 



These works underscore the relevance and timeliness 
of a better understanding of the dynamics of nonlinear 
waves in double well potentials. In view of that, in the 
present work we offer a systematic methodology, based 
on a two-mode expansion, of how to tackle problems in 
double wells, as regards their stationary states and the bi- 
furcations (and ensuing instabilities) that arise in them. 
This way, considering both cases of attractive and re- 
pulsive nonlinearities, we illustrate the ways in which 
a symmetric double well potential is different from an 
asymmetric one. In particular, we demonstrate that, con- 
trary to the case of symmetric potentials where symmetry 
breaking follows a pitchfork bifurcation, in asymmetric 
double wells the bifurcation is of the saddle-node type. 

The paper is structured as follows: in Section II, we 
present the model and set the analytical framework. In 
Section III, we illustrate the value of the method by 
highlighting the significant differences of symmetric and 
asymmetric double wells. Finally, in Section IV, we sum- 
marize our findings and discuss future directions. 

Model and Analytical Approach. In a quasi-lD setting, 
the evolution of the mean-field wavefunction of a BEG 
Q (or the envelope of an optical beam 5]) is described 
by the following normalized NLS (GP) equation. 
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In the BEG (optics) context, /i denotes the chemical po- 
tential (propagation constant) and s = ±1 for attractive 
or repulsive interatomic interactions (focusing or defocus- 
ing Kerr nonlinearity) respectively; below, for simplicity, 
we will adopt the terms attractive and repulsive nonlin- 
earity for s = ±1 respectively. Finally, in Eq. ^, V{x) 
is the double well potential, which is assumed to be com- 
posed by a parabolic trap (of strength fl) and a sech - 
shaped barrier (of strength Vq, width w and location xq); 
in particular, V{x) is of the form: 



V{x) = -n^x^ + V^osech^ 



X — Xq 

w 



(2) 



with the choice xq = {xq ^ 0) corresponding to a sym- 
metric (asymmetric) double well. Note that such a dou- 
ble well can be implemented in BEG experiments upon. 



e.g., combining a magnetic trap with a sharply focused, 
blue-detuned laser beam |2Q|. Similar double wells can 
also be implemented e.g., in optical systems. 

The spectrum of the underlying linear Schrodinger 
equation (s — 0) consists of a ground state, uq(x), and 
excited states, ui(x) (/ > 1). In the nonlinear problem, 
using a Galerkin-type approach, we expand u(x^ t) as. 



u{x, t) — co{t)uo{x) + ci{t)ui{x) 



(3) 



and truncate the expansion, keeping solely the first two 
modes; here Co,i(t) are unknown time-dependent complex 
prefactors. It is worth noticing that such an approxima- 
tion (involving the truncation of higher order modes and 
the spatio-temporal factorization of the wavefunction) , is 
expected to be quite useful for a weakly nonlinear anal- 
ysis. In fact, as will be seen below, we will be able to 
identify the nonlinear states that stem from the linear 
ones, as well as their bifurcations. 

Substituting Eq. Q into Eq. Q, and projecting the 
result to the corresponding eigenmodes, we obtain the 
following ordinary differential equations (ODEs): 

ico = {ujo - ^)ca- sAo\co\'^co- sB{2\ci\'^co + clco), 

- sri|ci|2ci-sro(2|co|'ci+cgci) (4) 
ici = (tJi -^)ci - s^i|cipci - sB(2|copci -t-CoCi), 

- sEolcopco - sri(2|cipco -t-CiCo). (5) 

In Eqs. Q-©, dots denote time derivatives, overbars 
denote complex conjugates, wo.i are the eigenvalues cor- 
responding to the eigenstates mq.Ij while Aq = J u^dx, 
A\ — J ufdx, B = J uffuldx, Fq = J uoUidx and 
Fi = J uiUgdx are constants. Recall that uq and ui 
are real (due to the Hermitian nature of the underlying 
linear Schrodinger problem) and are also orthonormal. 
Notice also that in the symmetric case {xq = 0), due to 
the parity of the eigenfunctions, Fq = Fi = 0. 

We now use amplitude-phase (action-angle) variables, 
Cj — Pje'^'^\ j = 0, 1 {pj and (j)j are assumed to be real), 
to derive from the ODEs Q-© a set of four equations. 
Introducing the function ip = (pi — (j)o, we find that the 
equations for po and 4>o are, 

Po = s[{ropl + Tipl)sm{ip) ~ plposm{2(p)], (6) 

4>o = {ujQ - p) + sAqpI + 2sBpl + sBpl cos(2(/7) 

+ s(J^ + 3popiTojcos{^), (7) 

while the equations for pi , (pi are found by interchanging 
indices 1 and in the above equations. Next, taking into 
regard the conservation of the total norm, we obtain the 
equation Pq + pf = N, where N = J \u\'^dx is the integral 
of motion of Eq. Q (the number of particles in BECs, 
or the power in optics). Finally, subtracting Eq. |(7J) for 
(po, and the corresponding one for (pi, we obtain: 

if == -Auj + s{AopI-AipI) 

cos{(p) 



sS(2 + cos(2(^))(p^-pf) 
[ropUpl-3pl)+ripl{3pl -2 



Equations ©, IJSJ is a dynamical system, which, in prin- 
ciple, can be thoroughly investigated using phase-space 
analysis (such an approach has been presented in [3,|3 for 
similar systems that were derived using different expan- 
sion of the field u). Here, we will focus on the fixed points 
of the system [corresponding to the nonlinear eigenstates 
of Eq. JQ)], and analyze their stability and bifurcations. 

Results. Below we will analyze all possible cases (s = 
±1, xq = 0, xo 7^ 0) for the double well of Eq. ^ with 
Vq = I, fl = 0.1, and w ~ 0.5 (the results do not change 
qualitatively using different values). 

First we consider the case of attractive nonlinearity, 
i.e., s = — 1, and a symmetric double well potential with 
xq = (implying that Fq = Fi =0). In this case, the 
parameters involved in Eqs. © and (jSJ are found to be 
Ao = 0.09078, Ai = 0.09502, B = 0.08964, wq = 0.13282 
and LUi — 0.15571. Then, it is readily observed that the 
possible real solutions of Eq. © are po — and pi = 0, 
as well as (/5 = (mod tt). The former two are contin- 
uations of the linear solutions in the nonlinear regime. 
However, the latter one is a non-trivial combination of 
the two modes for ip = n that results in an asymmet- 
ric pair of mirror-symmetric solutions |l5l Il7l | , emerging 
through a pitchfork bifurcation. From Eq. ((Sj), we ob- 
tain that this new branch of solutions bifurcates from 
the symmetric branch (po,Pi) = (VjV, 0) for 



N>Nr = 



At. 



35 - Ao ' 



(9) 



PoPi 

pI)]- 



(8) 



and for p < pc = ujq — AqN = 0.12115. These analytical 
predictions are in excellent agreement with the numeri- 
cal results Pc = 0.122(±0.001). It is also easy to see that 
the anti-symmetric branch {po,pi) = {0,^/N) does not 
give rise to such a bifurcation. The different branches 
of the full numerical solutions (including the bifurcating 
ones) have been obtained through numerical fixed point 
algorithms solving the steady state version of Eq. Q, 
and using continuation of the solutions over the param- 
eter p. The results are shown in the top left panel of 
Fig.l, where the norm of the solutions N = J \u\'^dx is 
shown as a function of the chemical potential (or prop- 
agation constant in optics) p. In addition, as expected 
from the nature of the bifurcation, the linear stability 
analysis has been used to illustrate the following: The 
emerging new asymmetric (i.e., "symmetry breaking") 
branch of solutions is stable, while the original symmet- 
ric branch is unstable beyond the bifurcation point due 
to a real eigenvalue A,, (see bottom left panel of Fig. 1). 
Next, in the same case (s = —1), we consider a dou- 
ble well potenial with a weak asymmetry (xq = 0.5). 
In this case, the constants involved in Eqs. lO-lO are 
found to be Ao = 0.14903, Ai = 0.15618, B = 0.02958, 
loq = 0.1249 and oji ^ 0.16535, while Fq = 0.0407 and 
Fi = —0.04077. Note that even such a weakly asym- 
metric case, renders the right well "shallower", in the 
following sense: the density, or power N (regarding the 
ground state of the linear problem), in the right well is 
smaller than the one in the left well. Thus, in the nonlin- 
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FIG. 1: (Color online) The top panels show the norm of the 
solutions of Eq.Q for attractive nonlinearity (s = — 1) as a 
function of n for symmetric (left panel, xq — 0) and asymmet- 
ric double well (right panel, xq = 0.5). The potential param- 
eters are Q = 0.1, Vb = 1 and w = 0.5. The solid lines denote 
the symmetric solution, the dashed-dotted lines denote the 
antisymmetric one, while the dashed lines denote the asym- 
metric solutions that are generated from the bifurcation at 
fj,c ~ 0.122 (pitchfork) and fic ~ 0.009 (saddle-node) respec- 
tively. The bottom panels show the maximal real eigenvalue 
associated with the linear stability of the symmetric branches. 



ear problem, the respective branches that bear the larger 
part of the density in the right or in the left well (i.e., 
the ones having, roughly speaking, the shape of a single 
pulse in each of the wells) are no longer equivalent. This 
results in a significant difference between the asymmetric 
and the symmetric case discussed above, namely there is 
no longer a pitchfork bifurcation, but instead, there is a 
saddle-node bifurcation. This result is shown in the top 
right panel of Fig. 1, where N is shown as a function 
of /x. It is readily observed that, due to the nature of 
the saddle-node bifurcation, two branches (one of which 
is stable and the other one is unstable) "collide" at some 
critical value oi ^ — ^c (see below) and disappear. These 
branches are the more "symmetric" one, that has sup- 
port in both wells (see solid line in top left panel of Fig. 
2), and the one pertaining to the state having the form 
of a single pulse in the shallower well (see dashed line in 
the rightmost top panel of Fig. 2). The instability of the 
former branch is depicted in the bottom right panel of 
Fig. 1, where the maximal real eigenvalue is shown as a 
function of /i. On the other hand, there exists also an- 
other single-pulse branch supported over the deeper well 
(see top third panel of Fig. 2), which persists all the way 
to the linear limit. Furthermore, the dash-dotted anti- 
symmetric branch of the top right of Fig. 1 is shown in 
the second top panel of Fig. 2. 

The novel feature described above, namely the asym- 
metric breakdown of the pitchfork bifurcation into a 
saddle-node one, is a particular feature of asymmetric 
double well potentials that, to the best of our knowl- 
edge, has not been previously appreciated. Notice that 
Eq. © predicts that the saddle-node bifurcation oc- 



FIG. 2: (Color online) The steady state solutions of Eq.Q, 
(see also Fig. 01 for the focusing, asymmetric case (top panels) 
and their linear stability (bottom panels) for /i = 0.05. The 
black-dashed line shows the double well potential. 
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FIG. 3: (Color online) Same as in Fig^but for the repulsive 
nonlinearity (s = +1). The dashed-dotted lines denote the 
symmetric solution, the solid the antisymmetric, while the 
dashed lines denote the asymmetric solutions generated from 
the bifurcation at fic ~ 0.168 (pitchfork) and /ic ~ 0.207 
(saddle-node) respectively. Contrary to the case s = —1, the 
bifurcations originate from the anti-symmetric branch. 



curs at fic = 0.08748, while the numerical result is 
He = 0.09 ± 0.001; apparently the two results are again 
in excellent agreement. 

Let us now consider the repulsive nonlinearity (s = 
+1). In this case, for the symmetric potential (xq = 
0), the pitchfork bifurcation still occurs; however, now it 
does not originate from the symmetric branch, but rather 
from the anti-symmetric one with (po,Pi) = (0, viV), 
giving again rise to symmetry breaking. Analyzing Eq. 
©, we find that this occurs when 



N > Ne> 



Alu 



3B-Ai' 



(10) 



and for jj, — ujo + 3BN = 0.16822, once again in excellent 

agreement with the numerical result /i^ = 0.168(±0.001). 

On the other hand, in the same case (s = +1) but 

for an asymmetric {xq = 0.5) double well, the bifurca- 
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FIG. 4: (Color online) The steady state solutions of Eq.Q, 
(see also Fig. |^ for s = -1-1 in the asymmetric case (top 
panels) and their linear stability (bottom panels) for /i = 0.22. 
The black-dashed line shows the double well potential. 



tion still originates from anti-phase (between the weUs) 
solutions, but again (as in the asymmetric case with 
s = —1), the bifurcation is of the saddle-node type. This 
is theoretically predicted to occur at ^c = 0.21342, once 
again in very close agreement to the numerical result 
He = 0.207 ± 0.001. The details of the bifurcation di- 
agrams, are illustrated in Fig. O (analogously to Fig. 1), 
while the steady state solutions and their linear stability 
are shown in Fig. 0] (analogously to Fig. 2). 

Conclusions. In conclusion, we have presented a sys- 
tematic analysis based on a Galerkin, two-mode trunca- 
tion of the stationary states of symmetric and asymmet- 
ric double well potentials. The analysis has been carried 
out both for repulsive and attractive nonlinearities and, 
as such, can be relevant to a variety of physical contexts; 
these include matter- wave physics (most directly), non- 



linear optics, as well as other contexts where it is rele- 
vant to consider double well potentials in the NLS model 
proper. We have demonstrated that our analytical ap- 
proach describes quite accurately, both qualitatively and 
quantitatively the features of the nonlinear solutions; nu- 
merical results were shown to be in excellent agreement 
with the analytical predictions. 

In the case of a symmetric double well potential, it 
has been shown that a symmetry-breaking (pitchfork) 
bifurcation of the ground state occurs for attractive non- 
linearities, while it is absent for repulsive nonlinearities. 
It has also been found that a similar bifurcation of the 
first excited state occurs in the relevant branches for re- 
pulsive nonlinearities, oppositely to the case of attractive 
ones, where such bifurcation does not happen. Addition- 
ally, regarding the above feature, we have illustrated that 
symmetric potentials are very particular (degenerate) 
due to their special characteristic of mirror-equivalence of 
the emerging symmetry-breaking states. We have shown 
that even weak asymmetries lift this degeneracy and lead 
to saddle-node bifurcations instead of pitchfork ones that 
were similarly quantified in both attractive and repulsive 
nonlinear! ty contexts. 

These results underscore the relevance of analyzing 
steady state features of nonlinear models (in the presence 
of external potentials) based on the states of the under- 
lying linear equations. It would be particularly interest- 
ing to examine the extent to which dynamical features 
of such models can be captured by similar truncations. 
Such studies are currently in progress. 
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